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ABSTRACT 


This thesis is concerned with proving that the eigen- 
moeules Ol a Specitic Unsymmetric matrix are real and positive, 
Peas nOuUceactUaliy Computing them. The method of finite 
PPh <enecoe to ool lIcCd Goeune vibration analysis of a 
cantilever beam and leads to an unsymmetric stiffness matrix 
Pmt Ceivelwmc Droplem formulation. The technique 
pimuovoosi tie Prool 1s based on @ perturbation theory 
Bien by VWalkinson for real symmetric matrices. Applicaticn 
of the theory is made to the cantilever beam eigenvalue 
Seeoleie tne results verity that the eigenvalues of this 
oOo vicr Vnaymnetric matrices can be proven real and 
BOsttive Wiuhout their actual values being Caiculaved. 

Piet variable Lormudation of the cantilever 
beam vibration analysis is also examined to illustrate that 
Duoirelhulty CONvainis all of the properties of a symmetric 


Mente 1c 
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ie ROU TAO 


The vibration analysis of many structures can be 
eceomplished by the application of finite differences to 
the partial differential equations governing the free 
undamped vibration of the structure. Assuming the motion 
Seebie SGructure CO be harmonic, with natural frequency @, 


Meeos Go the matrix equation for the eigenvalue problem 


i= wr PD (1-1) 


moere K ang M are square matrices and are referred to herein 
Pemene shleitiess mabvrax and the mass matrix respectively. 
The vontar A defines the natural made, When K and M are 
Sime trLrcswana K or M aS posturve devinite, the eigenvalues, 
a pacmuceteewoen born K and M are positive definite, 
ene is always positive [Rer. 1] . When K and M are developed 
temas enerey techniques, they are always symmetric; however, 
when they are developed for the same problem by applying 
ivemlmrereniCes tO the poverning partial differential 
equations of motion and the boundary conditions, they 
EredUcntly emerge numerically as unsymmetric, ealling ance 
question whether numerical results for age Wali wea Tact. be 
peal, 

In this thesis an UNSyMMNeLCLe Sta tiness Matrix. us 
developec by Mir eeLerereices Or une Cantilever beam. 


A perturbation theory for real symmetric matrices given 





by Wilkinson | Ref. 2] is used to establish that all of 

the eigenvalues of the unsymmetric matrix formulation are 
real and positive; hence the results obtained for the 
Poouremeiredqueneies are ~eal, representing undamped harmonic 


motion as assumed, 


10 





IT. DERIVATION OF THE MATRIX EQUATIONS 


Mae structure selected for analysis is the cantilever 


beam of length L shown in Figure 1. 


YA 


es ee uae es a 
I A Oc 
| ——-—— : 


BEAM CONFIGURATION 
FIGURE | 


Miemeauatron Soverning the free undamped vibration of the 


beam can be given in the form 


~ 


aio. NS 
ee . \ +2 


Pee wire ots tame, and 





(2-1) 


P= Y/L 
i = Ab 
[L = p/(Er*) 


where ¥ as the deflection of the beam, X is the axial 
coordinate , pts the mass density, EH is the elastic modulus, 


and r is the radius of gyration of the cross-section, 


Assuming harmonic motion, 


cp = cp sin wt (2-2) 


do 





Substituting Equation 2-2 into Equation 2-1 leads to 


mM 
a =\ @=5) 
in which A=[U « 


The boundary conditions on the beam are 


_ ag _ 
P= s = ai X =O (2-Ha ) 
and 
ach _ eee 


Let the beam be divided into ten equal intervais of 
feos Nias shown in Figure 1, The nodes j=-1, j=ll, and 
gale are fictitious nodes located off of the beam. Let 
the derivatives in Equation 2-3 and Equations 2-4 be 


BeerOxImMatec DV 


) 
eb. . 
‘i 


d 2 : 2 _ 
(S*) = { P 5-4 ey ae Hea) (2a) 


II 


(-b, 4 th) 7 (20) (2-58) 


Feast 
oF 
ie 

ee 
II 


: ee ae / (2h) (245 
( DP. ot Psy Pt Pag? (2h) (2-5c) 


1S} 
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> 


J 
mM 
d 

) 


LHe Vwileh i—-O te Teor eruncaulon error Of cach oL the finsce 


(D -h h- +h fn’ (2-54) 
ee Jen J Jae je 


O 
difference approximations in Equations 2-5 is O(h ). 


Oia Oia OMe = ea) y2, 5, ci. \ se. 6) 29 oo, LO, BQUaTIONS 


ili 





Peace | =Owweano Fquations 2-4b at j=10; and eliminating 
the variables ®4.Po->P1 > andPyy eventually leads to 


the matrix equation 


(eee OO . . wee « O P, D, 
“4“~ 6-4 10... . #0 P, P, 
Mm -+ 6-4 10... 0 |e, , 
O 1-4 6-4 10... 0 A, PD 
WeOele-P7O=4 1 0. O \ dein A (2-6) 
Oo . 01-4 6-4 1 0 0 Ly SP 
Oem 2 © T-4 6-4 1 O Pr Sr, 
O . .. O14 6-4 11 |B PR 
wes. OL HE 5 +2 | Le Ze 
oes. flO CH 2 | \Fi0 Po 


oo 


(ices No enc 2CentLty matrix. Note that the stifiness 
matrix in Equation 2-6 is not symmetric. However, if the 
last row of Equation 2-6 is divided by two, the resulting 
Momiiileatioy 16 Cymnetric. This results in the mass matrix 
Pem@omleieet Beane CCual to 1, Sinee the diagonal element 
teadow lace Ow 15 equal to one-half; however, M 1s st211 
eyumeurle and positave definite, and all cigenvalues are 
real. The symmetric stiffness matrix will be referred to 
berealtvor as K i Note that eal lof Ghe finite dilierence 


ite 


equations used to construct K. are centered about j. 


ill 


Lit can be shown that this symmetric formulation 
Homllcniurcain oO a Maprax equation developed using energy 
techniques with one-half of a mass interval at j=l10. 


1S. 





An alternate matrix formulation of the governing 
difference equations can be constructed in a similar manner 


by using the O(n=) approximation Ref. 3 





(S), Bee Py-07" Psy 208 13D, (2°) 
(250) 

fioswead of Equation 2.5c. 

omoarecenscecuence Of the fact that Hquation 2-/ is not 
eeeredsanOul jj, Che fictitious node j=l2 is not required, 
preigMeetom 2-5 1S applied ab j=1,2,35.<esee+-9. The 
boundary conditions are satisfied by applying Equations 2-4a 
Mee -O ond Equations 2-4b at j=10. Eliminating Pi; Pr > 


Pro; and PD, , leads to the unsymmetric matrix formulation 


mel Cid Dr ie 
| ao O., : : O PD, D, 
ieee or 6} 6lUUUC; O F, P| 
mn 
Oe ee tb + 6 4+ 1 0 0 OS ecu D. (2-8) 
O : OQ Ss eS are) as 2. 
O : ‘ Oe = eG 4 a Se, , 
O . . : O uli -3% 4s 17% % PD 
On . Oo #¢ -1 & Dy = 


The stiffness matrix in Equation 2-8 will be referred to 


as K ,. 
E 
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III. ANALYSIS AND RESULTS 


Peer cturbaLl1on theory given by Walkinson [Ref. 2 | InO1e 
fee cymme pric Matrices will be used to prove that the 
eigenvalues of the real unsymmetric formulation given by 
Equation 2-8 are real. The brief description of the theory 
Woasen TOllLows 1S expanded in the Appendix. 


Let 
= A + B 
Ko 


where he is the real unsymmetric square matrix of order n, 
Perse oneh xn real symmetric matrix, and Bis ann xn real 
Mie imecrie Mavrix.,. Every eigenvalue of Ky jies within at 
Heaeu one Of the n circular dises whose centers are the 


=] 


ry 


> e at ea Ee : peers , Hay 055) ee Pee: ac ay ae 4 a 
fee OV eS Nesb OL UM. ia oie eS ieee 6S oe ta Ce 


5 9, the e-norm of B, where 


ifell > = (maximum eigenvalue of BB) 


N]e 


and ef tenbaewiermitian COnjusate matrix of B. The matrix 
Peers wncm Hermmahian, the Cigenvalues of which are all 
Piel teaiyem Of the n circular discs form a connected 
COeiiudaGlmvead trom Che Overs, there are precisely m 
eigenvalues in this connected domain. In the case of real 
unsymmetric matrices, any complex eigenvalues must occur 
gewconjUeaue pairs. If the aoe CAreulLae dice 1Se cota ved 
and contains only one eigenvalue, this eigenvalue must 
Cia wore De reat. 

The matrix K. can be decomposed in many ways to obtain 


2 
the matrices A and B. One technique is to use the symmetric 


i 





matrix RY dees and obtain B by subtraction. This method 
was tried and resulted in such a large al 2 Baan Only 
two eigenvalues of Ky Cause Wem. Oem) oa  CoMsegquenuly. 
two other approaches were considered in an attempt to 
mae7e1Op 2 Smaller norm size that would prove all the eigen- 
values to be real and positive. These two approaches are 
Peoterred GO as Method I and Method II. 

ie ietned 1, theeiatrix K. 1S) decomposed such thag “5 


SeGerciiicm une smallest numoer of elements possible, i.e., 


eee O Ap -%3 Ls -1 


eo e O a =i 3, 
K 
Ee 
sal 6 aH lh O ® e e e e @ 
17-4 6-4 1 SarsOn 20.5 LO; 20 
-{- 
ia eR Po Oe 10 | 
| 
e Oy a) 1 -i ‘. e O -4 O O | 
a i 


BRquation 3-1 shows only that portion of SS that requires 
alteration. The eigenvaiues of a and the e-norm of B, 
were determined using the IBM SSP subroutine EIGENP and 
glemCapvietecra Tapile 1 and plotted in figure 2. Nove 

Hoven eimai vero tie Circular discs are isolaved 


eng a tne posavave domain. AS a result, applying the 


perturbation theory proves that at least five of the nine 


16 





eigenvalues are real and positive. Also note that the 
Pove real Cigenvyalues are the five largest eigenvalues. In 
order to prove that the lower eigenvalues in the overlapping 


Meacwlareciscs are real, the inverse formulation 
1 i haves Oe 
is examined. If “aN can be proven real and positive, then 


\ is real and positive. 


IDNs ee pal 


KIGENVALUES AND NORMS - METHOD I 


Eigenvalue ): 


t+ 
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3 
a 
> 
6 
¢ 
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, 





Value of ie =0.814 


r|l2 


Value of |B =180.9 


| 2 
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e S3yndis 
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vA 
608! = {ial 


“y 4O SSNIVANZOIZ X 
ty 430 S3NIVANZOIS | 


3u 


008 OOL JONES 


O0& 002d 


OO! 


ANWId-X3 IdWOD SHL Ni SLINSSY I GOHLAN 
ye | sl DRO ee 





PwieOomilveroion darectly relates the lowest eigenvalue of 
ale 
Ky to the largest eigenvalue of = Pe ond (SO. On. oAnce the 


meet —elgenvalues of K,. have not yet been proven real, the 


& 


mesier e1senvalues of KS are of particular interest. 


ihe decomposition of KT ieee Wkeeele\vorel ya iiea bs: 


nN 
= 
Nth 
Re) 
iD 
NY) 
OD 
Ce) 
Ne 
te 
= 


LI Ls, 2s 
Soe eo 
ee S205 s5 46. S54 67 


fem 
Co 
LO WN 
Ne 
eI 
—- N 
od 
es) 8, 
NY) 


LO N WW 


2 

p ee Ome.) 55 Bers 60 145 c= 
Beiyy23 368 558 73 9h 109 1992 

| 34 13 275 46 675 91 1155 140 259 

y ie ce ee cO)6OS eer ay? 322 

(4b 17 365 


62 925 127 lols eon 367 | 
cae 
2 
=) a io 4S 
lL ; me ies: - 217 
32 363 
54 662 
a 
80 924 
109 127 | 
a eer. 264 


4 15 32 54 80 109 140 172 204 
45.17 363 62 925 127 1641204 387 





Pt 
He 


20 





: (3-3) 
| 
| 


= =) © © © © (2) © © 
© @ © © ®) S © = 
Oh 
NO 
Nf} 


ee 


ioewcifpenvalues or Aa and the @=-norm of B , are also 


cae@ywiotoeteable £ and in Figure 3. Examination of Figure 3 


pa | 
reveals that the iargest elgenvalue of oS is real and 


Pee mervesmnemee, the lowest eigenvalue of K, is real and 


a 
Beectumve. lnls, Six Of the nine eigenvalues have been 
ELOvemerecaleand positive. To verify that these eigenvalues 
emer rcal amd positive, the actual eigenvalues of 

K, and — were computed using EIGEMNP and are tabulated 
iijmiaote mt amd pHLOLLed ifn Figure 2 and Figure 5. Iney do 
Mmio ker Witiin the circular discs. 


The decomposition of Kp and KS 


Or Mepmod jl led vora 
Cowive Wpellorse morm and the proot Ghat only six of the 
(ite memreMyalves were real and positive. Consider now 
Mecnod = Lila which K,, is decomposed such that the moduli 
of the elements of B are as small as possible (making 


ali 





B skew-symmetric), i.e., 


(3-4) 
6 aL 4 O 
ey! ons 1 = 
2 -34, 4s -1>/1, 
5 ge 2 ol 
> 
: eet il O O C0 O 
ie 3 
ee 6 “3 Ji O O ~+79 t 
S: 
0 34 Jy | / 0 -b Fg 
Alt one 


wuyeveMaeemoOrtion Of K. requiring alteration is shown in 


Zz 
Equation 3-4. Using the same procedures as in Method I 


Meco tO tse results tabulated in Table IL and plotted in 
Figure 4, The norm of Bie is much smaller than the norm 


of B_, and examination of Figure 4 reveals that seven of the 


ee 
Dime weilrenvalues Of K aconleais aiumpOcsuiwe,9 DEConnos ime 


: in a manner similar to Equation 3-4 leads to a 


and Bly! The results are also shown in Table II and in 


Figure 5. Using this formulation, the lowest eigenvalue 


K 


of K, is proven real and positive. Consequently, use of 


2 
Method II proves eight of the nine eigenvalues to be real 


Ce 





and positive. The last eigenvalue must also be real since 

ieecam have no Conjugate. It must also be positive since it 
is larger than the lowest eigenvalue which has been proven 

Peglulve, wthereitore all the eigenvalues of Ky are real and 
os itive, 

A third method considered is to treat K, in a manner 
Similar to the conversion of Equation 2-6 to the symmetric 
formulation Ke ide wa lSeime piaod,. K, ie sole 6 red aby ml tees) lyr ie 
row 9 by two. Although this leads to another unsymmetric 
mabrix, 4 much smaller |BI| 2 Cam be developed.” However, 
from past experience it is concluded that this approach 
feulaenet prove all eigenvalues real and positive without 
hecoruane Co the inversion of the matrix. Since Method I1 
Beers he Picenyalues of QO real and positive without 
miPerautOn Ol Lae matrix, and since the eigenvalues of This 
third matrix formulation are not Lt this approach is not 
pursued Peer , 

As a check on the accuracy of the calculated eigenvalues 


of K pmoonparicon Onl The Cxace natural Trequencies | Ref. | 


2 
and the numerical natural frequencies is made in Table IIil 


Pemeuee ave lowest frequencies. 
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TABLE IT 


BKIGENVALUES AND NORMS - METHOD IT 


Eigenvalue ): 


thos Ke Aves Ken 
1 -O.144 O0013 -5 58 0.066 
2 0.040 0.049 Ono7@ 0.079 
3 0.363 0.369 0.0844 Oy aL OF 
ly alge oF on 
> Se O. On 
6 oe O. OF 





Value of rz ||. = 0.148 


Value of Par || a 90.3 


a4, 
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Sbr1'o = {lal 
fy 40 SSNIVAN29I9X 
Dy 40 S3aNIVANSDIS } 
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TABLE III 
NATURAL FREQUENCIES OF CANTILEVER BEAM 


Natural Frequency rad/sec 


Exact Ref. ul Numerical 





at, 





IV. MULTIPLE VARIABLE FORMULATION 


An alternate unsymmetric formulation is obtained when 
the two variables @ and ™N where /) is the bending moment 
in the beam, are used instead of only D. The governing 


equations are 


=P _ Ip (4-1) 


\ 
x= 





Cr 


and 





Le. =p (1-2 





mew moundary Conditions for the two-variable formulation 


elias 
dp 
Cee) O 
gcse é ea 23 a 
W\wu J EX Ve Si aa Me 
and 
dy = 2% 
= ——- =0 


Consider the beam shown in Figure 1. Applying Equation 
Paimomomucuacionm 4-2 at j=1,2,35 e0e.sese- 10, Equation 4-3 
at j=0, Equation 4-4 at j=10, and eliminating D>P.> 


P22 MNo: and? 14 leads to the matrix equation 


28 





(=>) 


a2) -né 


O 














Nae 











ase, 





Equation 4-5 is unsymmetric. However, if the variables 27,, 
Yo: 5 0h 3s ee ee 7h 10 are eliminated, Equation 4-5 reduces 
to Equation 2-6. As previously discussed, Equation 2-6 is 
Simply a symmetric matrix in disguise, so that the eigenvalues 
are real and positive. 

The two-variable formulation has been presented here to 
point out that many unsymmetric matrices may in fact be 
Recneibie FO Symmetric matrices; hence, they have all of 


the properties of symmetric matrices. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


Micom@astemeconccpe Litistratvead a0 this thesis is that 
mone f2iven UnSyMMeurile Matrix the elgenvalues may be 
emorined for proot that they are real and positive without 
Peuwuelty calcuiating them. Use of the perturbation theory 
given by Wilkinson gives a relatively straightforward method 
by which this may be accomplished. For the simple example 
petleevea Caleculavion of the eigenvalues is easily carried out, 
eve CONltirming the results. 

(ieee conmuecnidea thay this technique be applied to other 
elecmvealve problem formulations of more involved physical 
Svetems. If the technique proves itself reliable in appli- 
Salou coenabraces Of large order, significant computer time 
can be saved by ascertaining whether a given modelling method 
will GCootMmiieGecl, DOsitive elgenvyalues prior to conducting 


Hen oc eeisenvalue computations. 


Bik 





APPENDIX 


PERTURBATION THEORY GIVEN BY WILKINSON 


With slight changes of notation and text, the following 
Pitteticielemerecaroqucea for case of reference from the pages 
Pated for Reference 2. 

Consider an eigenvalue A (possibly complex) of the real 


_ (A + B)x =.Ax (A-1a) 


fHeewiech A is symmetric but Bis not. Since A is real 


Biere 1S an orthogonal matrix H such that 


= dias r, ) (A-1b) 


and the A. are real. Because A is an eigenvalue of (A+B) 


— 


the matrix (A+B-AI) is singular and hence its determinant 


is zero. From Equation A-la and Equation A-1b 


HO (ALB AI)H = diag O47 A \4+H BH (A-2) 


pWoOstermmeowcdeterminants, the matrix on the right of Equation 


Pee Mcteealso be Singular, Two cases are distinguished. 


Case 1. A=), for some i 


Case 2. NARA for any 1, So thay 


=k 


diag(),- \ 4H BH = diag( rar ) | Tedtag( ae) H BH | (A-3) 


and taking determinants again, the matrix in brackets must 


be singuler. If (I +Z ), whereZis a general matrix, is 





Sno eles, \z | 9 21; Tor if |[z o<; none of the 


Se 





eigenvalues (I + Z ) can be zero Ree. 5 | me ecnce 


Jaieg (Q,- \ 74H BH lL. AG (A-L ) 
Sua Vey es 


mec DA NLT Helle lle 22 os) 


i.e. (in both Case 1 and Case 2) 


min |\,- | 4 ll Le Pll. (A-6) 


mor at least one ), for every eigenvalue , . That is to 








say, every eigenvalue of (A + B) lies in at least one of the 


secular = dises siven by 


PAL s [Po 








cffellll. «7 





and since H is orthogonal (all eigenvalues of which have 


rnetgjawbas ~ <- a s4M¢ + 
Pee A ec, OO OHhe pis bhesp Ted ICA SD] 





eal =lll, 


Wwiw@enemit any S Ol vhnese discs form a connected 
Moiemimtooteued from the otners, then there are precisely 
Seomoenverues 2m this domain. 

In the case of a real unsymmetric matrix such as (A + B), 
any complex eigenvalues must occur in conjugate pairs. If 
hi-r,| Se | 


isolated and contains only one eigenvalue. This eigenvalue 


B 





o (J # i), then the iv” circular dise is 











ue veui he woOrempe roal. sHence, if all eigenvalues of A 


are separated by more than 2 





|B] |o: all eigenvalues of 


(A + B) are real. 
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